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Chapter 1 
Introduction 


1.1 General Introduction and Objectives 

Structural components used in aero-propulsion applications are often subjected to 
complex multiaxial stress states while at high temperatures. Under such demanding 
conditions materials can incur permanent deformations and changes in material state. 
When faced with such a difficult situation the engine designer would be greatly aided by 
a viscoplastic multiaxial deformation model that could accurately describe the material 
response over a wide range of conditions. Viscoplastic models such as the GVIPS 
developed by Arnold et al.[15] or Bodner-Partom [18] models can perform this task. 
Unfortunately, these models have yet to be validated under complex service conditions. 
This is due to the fact that experiments used for validation must subject the material to 
multiaxial stress states as well as elevated temperature. In short, the testing is very 
complex, expensive, and time consuming. 

The work presented here is an extensive study of the evolution of yield surfaces 
after axial-torsional prestraining. It consisted of two phases, the experimental research 
conducted at NASA Glenn Research Center (NASA-GRC), and data analysis conducted 
at Penn State University. The goals of the study were: (1) build a library of yield surface 
data for Inconel 718, Haynes 188, and 316 stainless steel, (2) mathematically fit the yield 
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surface data with a model incorporating yield surface distortion effects, and (3) use the 
model parameters to construct evolution equations describing the yield surfaces. 

1.2 Plasticity Theory 

Classical plasticity is a mathematical theory that describes permanent time- 
independent deformation of materials, such as metals at room temperature. Such 
mathematical theories are phenomenological in nature and are based on experimental 
observations. Classical plasticity theories have three distinct parts, a yield criterion, a 
plastic flow rule, and a hardening rule. The current section will discuss these three topics 
as well as loading criteria and will then introduce viscoplasticity. 

1.2.1 Yield Functions 

Yield functions are used to separate elastic stress states from those where plastic 
deformation occurs. For an isotropic metal under isothermal conditions the general form 
of the yield function can be written as, 

f = F(<j lj )-k = 0 (1.1) 

where F(Oij) is a function of the current stress state and k is a constant based on the initial 
tensile yield strength of the material. Yielding initiates when f=0 and the surface defined 
by this condition is know as the yield surface. All stress states inside of the yield surface 
are elastic and stress states outside of the yield surface are not permitted unless a 
viscoplastic model is being used. As a result, the yield surface must evolve as plastic 
deformation occurs in order to insure that the stress state remains on the yield surface. 
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The two most popular yielding criteria are associated with Tresca [1] and von 
Mises [2] criteria. Both of these theories are for isotropic materials and are deviatoric in 
nature. 

The Tresca, or maximum shear stress criterion is based on the assumption that a 
material will yield if the maximum shear stress exceeds the critical shear strength, T> . 

The critical shear strength is defined as Xy = a y /2 where g y is the tensile yield strength. 

In terms of the principal stresses the Tresca yield condition takes on the following form, 

/ =max(|cT I -a n \,\a n -<t ; |)-(t v (1.2) 

where Ci, O//, Ow represent the principal normal stresses and oy the tensile yield strength. 
Equation 1.2 assumes a fully three-dimensional stress state, however the work in this 
study involves only stress states in the axial-shear stress plane. Consequently, there are 
only two non-zero stress terms: axial stress, On, and shear stress. On- Therefore, the 
Tresca yield function reduces to, 

07, + 4<7, 2 2 = 07 (1*3) 

The Tresca yield criterion is very easy to apply, however it often provides a conservative 
estimate of the yield stress. 

The von Mises yield criterion often agrees better with experimental results than 
the Tresca criterion. This criterion is also known as the distortion energy criterion 
because it predicts yielding to occur when the maximum distortion energy exceeds the 
energy required to cause yielding under pure tension. The von Mises theory can be 
expressed as, 

=0 (1.4) 
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where J 2 is the second invariant of deviatoric stress given by. 


J 2 


= T S U S U 


(1.5) 


and deviatoic stress, sy, is defined as. 


S U G 'j <5 & kk^ij 


( 1 . 6 ) 


where 8,y is a second order identity tensor. When expressed in the axial-shear stress plane 
the von Mises yield criterion becomes 

<7n+3 a,2=<r*. (1.7) 

Figure 1 . 1 shows the relative size of yield surfaces predicted by the Tresca and von Mises 
yield criteria. Also note from Esq. 1 .7 that if the shear stresses are multiplied by V3 the 
von Mises ellipse becomes a circle. This characteristic will be used extensively in later 
analysis. 

1.2.2 Loading Criteria 

In order to discuss plastic flow it is first necessary to present a formal definition 
of loading. There are three different types of loading that can occur in a work hardening 
material when the stress state is located on the current yield surface. If an infinitesimal 
stress increment is added to the stress state directed outside of the yield surface then the 
loading condition becomes. 


/ =0 and 




d<*y> 0 - 


( 1 . 8 ) 


However, a stress state located outside of the yield surface is not permitted therefore the 
yield surface must evolve to accommodate the new stress state. 
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The second loading condition is one where the stress increment is directed such 
that it moves the current stress state inside of the yield surface. This is called unloading 
and is given by, 


/ = 0 and 




(1.9) 


The final loading type is called neutral loading, which occurs when the 
infinitesimal stress increment is directed tangent to the yield surface. This type of 
loading causes no evolution of the yield surface and can be represented by 

f =0 and — do 7 = 0. ( 1 . 10 ) 

da :i 1 


1.2.3 Flow Laws 

As already stated, when a material reaches a stress state where f=0 it ceases to 
behave in a linear elastic manner and incurs permanent deformations. Because of this 
permanent deformation the material is said to undergo plastic flow. This plastic flow 
changes the material state, which drives the evolution of the yield surface. Since plastic 
strain is path-dependent, it must be computed in an incremental form, d£jj. 

It is often convenient to define the flow law in terms of a plastic potential function 
such as the one presented by von Mises [3], 

de?=dX^- ( 1 . 11 ) 

where Q is a plastic potential function that is a scalar function of stress, dX is a 
nonnegative scalar that is zero unless the stress state is on the yield surface and the 

loading criteria is satisfied. It can be shown that ' s normal to the plastic potential 
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function, £2, therefore the plastic strain increment is always normal to the plastic potential 
function. Consequently, Equation 1.11 is known as the normality flow rule. 

If the potential function in Equation 1.11 is a yield function then the flow law 
becomes an associated flow law. If the von Mises yield criterion is used (Equation 1.4) 
then Equation 1.11 is known as the Prandtl-Reuss flow law (Prandtl [4] and Reuss [5]), 
given by, 

de? = dAs,j . (1.12) 

1.2.4 Hardening Laws 

Hardening occurs when the conditions in Equation ( 1 .8) are met. Since the stress 
increment tries to push the stress state to fall outside of the yield surface the surface must 
evolve. This evolution is called hardening. There are three types of hardening that are 
relevant to this study. They are isotropic, kinematic, and distortional hardening. 

Isotropic hardening occurs when the yield surface increases in size without a change in 
shape or translation of the center point. This type of hardening implies an overall 
increase in the yield strength of the material. In contrast, kinematic hardening is a change 
in the location of the center of the yield surface without a change in size or shape. The 
Bauschinger effect results from kinematic hardening. Figure 1.2 shows changes in the 
yield surface due to isotropic and kinematic hardening effects. Distortional hardening is 
more difficult to represent mathematically than the other two types. It is a change in the 
overall shape of the yield surface. This often involves a change in the curvature of the 
yield surface, usually an increase in the front and a decrease in the back. An example of 
curvature change after an axial prestrain is shown in Figure 1.3. Another distortional 
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effect is known as a cross effect which is a change in the size of the yield surface normal 
to the direction of prestraining. Positive cross effects are shown in Figure 1 .4. 

A general yield function, such as the one in Equation 1.1 can be easily altered to 
incorporate isotropic hardening effects. This is accomplished by making k a function of a 
state variable q that accounts for the material loading history. As the total plastic strain 
imparted on the material increases, k(q ) increases, resulting in an overall increase of the 
material yield strength in all loading directions. 

The yield function in Equation 1 . 1 can also be altered to represent kinematic 
hardening as follows, 

f = F{a, (1.13) 

where a ,y is known as the back stress tensor. The back stress describes the translation of 
the yield surface with respect to the plastic deformation history of the material. Note that 
k is a constant. 

It follows that these two hardening types can easily be incorporated into a single 
model to describe both isotropic and kinematic hardening as shown, 

/ = F(Ojj — GCy ) — k (q) (1.14) 

This mixed hardening model often better predicts yield surface evolution than either of 
the previous two, however, it still does not account for any change in the shape of the 
yield surface. 

1.2.5 Distortional Hardening 

Being a change in shape, distortional hardening is not easily described 
mathematically. The change in shape is greatly influenced by the direction of loading. In 
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order to discuss distortional hardening, a preloading direction must first be defined. The 
direction of the plastic loading is known as the preloading direction and dictates the 
nature of the distortional effects seen in the subsequent yield surface. Figure 1.3 shows 
yield surfaces with a tensile preloading direction. 

An important point is that the prestrain and prestress directions are not always 
coincident. This is particularly true for nonproportional loadings. For example, let a 
specimen be loaded in tension into the plastic regime under strain control. If loading is 
then continued, under strain control, in the shear direction the stress will not follow the 
same path as the strain. Once shear straining is initiated the axial stress will begin to 
reduce due to plastic coupling effects. However, the axial strain does not change when 
shear strain is applied. This is illustrated in Figure 1.5, where figure (a) shows the strain 
path followed while figure (b) show's the response of the stress. The arrows show the 
direction of the change in stress at the given point. 

There have been many approaches used to describe distortional hardening effects. 
One approach, by Ortiz and Popov [6], uses trigonometric functions. The yield function 
that they used is, 

f - | - a Y (l + p 2 cos 26 2 + p ? cos 30 j ) ( 1 . 15 ) 

where IL%II is the norm of the effective deviatoric stress defined as = ($y - a y), ay is the 
initial tensile yield strength, p 2 and p3 are model parameters that describe cross effects 
and changes in curvature respectively, and 0 2 and 63 give the angle between the effective 
deviatoric stress tensor and the phenomenological internal variables p 2 and P3. These 
internal variables are unit tensors that define the direction of the yield surface 
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characteristic being modeled by that term. For example, since the term with subscripts 3 
represent distortion, P 3 would represent the direction in which the distortion occurs. This 
model has the ability to represent isotropic, kinematic, and distortional hardening, 
however the model parameters have little physical meaning. 

Another model developed by Eisenberg and Yen [7] uses a second order tensor to 
describe distortion effects. The yield function for this model is, 

/.It, + <11‘> 

where fly is a tensor that describes distortion. The model was shown to accurately 
reproduce isotropic, kinematic, and distortional hardening effects. The parameter fly 
models both distortional and isotropic hardening of the yield surface. The rate equation 


for fl,y is. 


fl. 


= -Am , 7 


( 1 . 17 ) 


where wy is a unit tensor describing the direction of distortion and A describes the 
magnitude of distortion. Isotropic hardening is handled by defining fly as proportional to 
the initial yield stress. 

Zyczkowski and Kurtyka [8-10] constructed a model to describe yield surface 


distortion based on geometric arguments. The foundation of the model is based on the 
ability to describe a closed surface from two circles with different radii and center 
locations. The model is given by. 


a .i = + [fern'* f - d u)k d u)k + R u> 1 + a j 


( 1 . 18 ) 
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where Qij is a direction tensor, /?, is a vector containing the initial radii of the yield 
surface in the 5 independent dimensions, d k measures the amount of distortion in each 
dimension, and t k gives the direction of the distortion. It should first be noted that the 
subscripts /, j, and k, range from 1 to 5 due to the use of Ilyushin’s 5-dimensional stress 
space. This model also does not follow the normal format of using a yield function, 
instead it directly calculates stresses at yielding. This means it does not follow with 
normal assumptions made in plasticity theory. 

Rees [11] also developed a distortional model that uses up to rank six tensors to 
fully describe material hardening. One form of this model is given by, 

/ = Cue,; + C tj a k , + C ijklmii G :j a kl (T mn - k ( 1 . 19 ) 

where C tJ , Cyki, and are tensors whose parameters are characterized from 

experimental data. Consequently the model becomes very cumbersome to use. 

A much more manageable tensor-based distortion model was developed by 
Voyiadjis et al. [12-13]. This model uses a rank four distortion tensor following Hill [14] 
and requires only four parameters for characterization. The model is defined as follows, 

/ = M m ( s u ~ a ij X s » " )- k 2 (1.20) 

where. 




^ + ^y>ikSjl +S il S jk )+UyU W |c- //"’ + -O (1-21) 


where a, j is the deviatoric back stress defined as. 


1 s 

a ij a ij ^ ^kk djj 


( 1 . 22 ) 
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In the model, a is a parameter defining isotropic hardening, c and d define front and back 
distortion respectively. The tensor t)jj is a unit tensor providing the preloading direction, 
where a unit tensor is given by Il'Ujjll = 1. Furthermore, ft ) and ft represent Heaviside 


functions which are defined as, 


H'"=H(v„(s n -a„)) 
H (2) =H(-v„(s n -a n )) 

where the Heaviside function is defined as, 

[1 if jc > 0 


( 1 . 23 ) 


H(x) = 


0 if x<0 


( 1 . 24 ) 


This model has the capability to be linearized which makes it very useful to for this work. 
This feature allows the model parameters to be found using linear regression. This 
process will be discussed further in section 3.3. 


1.3 Overview of Viscoplasticity 

Time-dependent material response can occur in nearly all structural materials, but 
it is especially important at elevated temperatures. As a result, a model used to describe 
deformation in time-dependent material must be able to handle the effects of phenomena 
such as creep and stress relaxation. Models based on classical time-independent 
plasticity often have separate equations for time-independent and time-dependent strains. 
Making them seem like completely unrelated mechanisms. However, unified 
viscoplasticity models attempt to include all permanent deformations into a single 
inelastic strain term. The following sections briefly discuss the thermodynamic 
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framework that can be used to develop viscoplasticity models. Two viscoplasticity 
models are then summarized. 

1.3.1 Thermodynamic Basics for Unified Viscoplasticity 

Elastic and inelastic deformations in materials can also be thought of as reversible 
and irreversible processes. From this, the material response can be represented in a 
thermodynamic framework. When constructing a viscoplasticity model from 
thermodynamic arguments the basic principals must be obeyed. 

• Conservation of Mass 

pdV = constant (1.25) 

where p is density and dV is an infinitesimal volume element. 

• Conservation of Linear Momentum 

J* pb : dV + j r)R o^rijdS = f R pa t dV (1.26) 

where b are body forces, a are accelerations, dS is an infinitesimal area, and n, are 
the components of a unit vector normal to the surface. 

• Conservation of Angular Momentum 

l P e ijk x jb k dV + l^e^XjO^dS = \ R pe ijkXj a k dV (1.27) 

where xj is a position vector and e&\ is the permutation tensor. 

• Conservation of Energy 

pu=O ij £ ij +pr-q u (1.28) 

where it is the internal energy density, r is heat supplied, and q u heat flux. 

• Clausius-Duhem Inequality 
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PS-P- + 




V J J 


>0 


(1.29) 


where s is the specific entropy and T is the absolute temperature. 


1.3.2 State Law 

In viscoplasticity the material state is more complex to define than in rate- 
independent plasticity. It is common to not only define the stress state (Gy), but also the 
absolute temperature ( T ) and an array of internal variables (£«) to characterize the 
material state. Note, ^ can be a combination of vectors and tensors. The thermodynamic 
potential is given by the Gibbs free energy, 

G-C^e^-H (1.30) 


where H(Zy, T, is the Helmholtz free energy given by, 

H = u-Ts (1.31) 


where u is the specific internal energy and 5 is the specific entropy. The Gibbs free 
energy can be differentiated to give the following. 


dG 


dG dG dG 

dZ ii + dT K da 


(1.32) 


such that. 


dG 

eU ' dO,: 


dG 


s = -- 


Pa = - 


dG 


dT’ dC„ 

where p a is a generalized force term. From this the total strain rate is given by. 


(1.33) 


d 

£>i ~ dt 


' dG A 


V 


dCT, 


d 2 G ^ d'G + d'G £ 

it * — \ -nk ba 


do -do ikl k ‘ do tj dT dO'jdi; a 
Under isothermal and linear elastic conditions the second and third terms go to zero. 


(1.34) 
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1.3.3 Dissipation Potential 

With the material state defined it is now possible to derive a flow law for 
viscoplasticity. In rate-independent plasticity the flow law was written in terms of a 
plastic potential function. Here a similar potential function is used, called the dissipation 
potential, Q(0,j, T, £«). As a result the flow law can be written as. 


_ dQ 
~ d(T o ' 


(1.35) 


It is clear from this definition that the dissipation potential evolves as inelastic strain is 
accumulated. The evolution equation can be expressed in terms of the internal state 
variable such that. 


Pa 



(1.36) 


where p a is the first derivative of p a with respect to time. When isothermal conditions 
are considered, Equation(1.29) can be given in terms of the Gibbs free energy. 


Pa 


d 2 G 

KM, 



(1.37) 


1.3.4 GVIPS Unified Viscoplasticity Model 

The Generalized Viscoplasticity Model with Potential Structure (GVIPS), 
developed by Arnold et al. [15-17], provides a good example of a model derived from the 
method discussed in section 1 .3. 1 -1 .3.3. This model uses a yield criterion, one internal 
variable (back stress, a jj), and an evolution law to account for nonlinear hardening. 

For the GVIPS model the Gibbs free energy is given by. 
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( 1 . 39 ) 


G = G E +G,+G A =~C iJkl o ij a kl -o^ -B 0 {g + B lg p ) 


where e'"is the inelastic component of strain. The dissipation potential is. 




q+l 


n + 1 q + 1 

For the GVIPS model the internal state variable is the back stress, a ij, which has the 
conjugate Ay. The model is governed by three basic equations: the flow law, 

K = 


the evolution law. 


A 


Qi Jkl C klpq b pq ifa fl Z,<0 




if a I >0 

IJ IJ 


and the internal constitutive rate law. 


where 


(l ij — Ljj kl A kl 


A = . 


1 IJ IJ 


27, 


L ijk , = G,J =K l {i ijkl +K 2 a iJ a kl ) 

b, = ~ *3 


/ = ( — — y. 


K 




( 1 . 40 ) 

( 1 . 41 ) 

( 1 . 42 ) 

( 1 . 43 ) 

( 1 . 44 ) 

( 1 . 45 ) 

( 1 . 46 ) 

( 1 . 47 ) 

( 1 . 48 ) 
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In the above equations, /is the yield function, y defines the stress below which only 
elastic strain exists, 7,y*/ is the fourth order identity tensor, ly is the effective stress, 7 2 is 
the second invariant of deviatoric stress, and ( ) are the MacCauley brackets, defined as. 


JO if x < 0 
jjc if x >0 


(1.50) 


The constants K\, Ki, and K 3 contain the model parameters (Bo, By, R a , /co, k, fi , n, p , q). 

1.3.5 The Bodner Partom Model 

The Bodner-Partom model (Bodner and Partom [18-19]) is presented as an 
example of a viscoplasticity model that does not follow a thermodynamic framework. 
Instead the model is based on dislocation dynamics and no formal yield function is used. 
The Bodner-Partom model allows inelastic strains under any loading condition. 
However, there are many loading conditions where the inelastic strain is insignificant 
compared to the elastic term. A basic outline of the model is provided below, for more 
information see Bodner and Partom [18]. 

The total strain rate is first broken down into an elastic and inelastic part, 


£ = £ e + £ in 


(1.51) 


where e' and e” are the elastic and inelastic strain rates, respectively. The elastic strain 


rate is given by the time derivative of Hooke’s law and the inelastic strain rate is given by 
the Prandtl-Reuss flow law (Equation 1.12). In this case the constant X accounts for 
isotropic hardening and is given by, 
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( 1 . 52 ) 


A = 


Va 


exp 


J>o 

( Z~ \ ( n + 1 


v 37 2y 


In 


where D 0 , Z, and n are model parameters. Kinematic hardening is represented by using 
an effective internal variable. 


Z eff = Z 0 +q\Z{r)dT + {\- q\ J Z(r>v dx 


( 1 . 53 ) 


where 


r ij = 




are the current stress direction cosines and the evolution equation is given by. 


Z = 


m(Z, -Z)W in 


( 1 . 54 ) 


( 1 . 55 ) 


where W in = o tj £ is the inelastic power. The parameters that must be determined are 


Z 0 , Zi, D 0 , m, n, and q. 


1.4 Yield Surface Experiments 

There have been many experimental studies on the yield characteristics of metals. 
In this section, an attempt is made to review previous yield surface work that most 
closely pertains to the current study. Excellent review papers on yield surface 
experimental studies were written by Hecker [20] and Michno and Findley [21]. 
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1.4.1 The Definition of Yielding 

In most yield surface experimental studies the specimen of choice is the thin- 
walled tube. This type of specimen can easily be subjected to various ratios of axial and 
shear stress in order to map a yield surface in the axial-shear stress plane. In addition, the 
specimen can be subjected to internal pressure in order to define a three-dimensional 
yield surface. 

When mapping yield loci there are two basic approaches. First, use a separate test 
specimen for each point on the surface, or use the same specimen to map an entire yield 
surface. The second approach is the most common because of the cost benefits and the 
elimination of specimen-to-specimen scatter. However, when using a single specimen to 
determine multiple yield points the state of the material must not change between each 
probe. This leads to a careful consideration of the definition of yielding. 

There are many different definitions of yielding ranging from the standard 0.2% 
offset strain rule, which accumulates 2000 pe (p£ = 10' 6 m/m) of plastic deformation, to 
the proportional limit definition, where there is no plastic strain accumulated. The three 
most common definitions of yielding for yield surface evaluations are the proportional 
limit definition, the offset strain definition, and the back extrapolation definition. Each of 
these definitions can produce significantly different test results due to the sensitivity of 
yield to small amounts of plastic strain. Representations of these yielding criteria with 
respect to the stress-strain curve are shown in Figure 1 .6. 

The proportional limit criteria defines yield to occur at the point where plastic 
strain begins to accumulate. This method requires very precise measurements of strain in 
order to be sure that yielding is detected as the onset of non-lineanty in the stress-strain 
curve. Phillips et al. [22 - 25] used a version of the proportional limit definition in 


NASA/CR— 2001-211162 


18 



combination with a back extrapolation technique in order to define yielding with zero 
plastic strain. The method used involves first loading until two consecutive data points 
deviated to the same side of the linear elastic loading line. From here the last three points 
were fit with a line and yielding was defined as the intersection of this line and the linear 
elastic loading line. 

The offset strain method is considerably easier to implement in experimental 
investigations. However, the magnitude of the offset strain is largely arbitrary, but 
different values can be required depending on material and loading rate. When used in 
multiaxial experiments the offset strain is given by an equivalent offset strain, such as. 



where ef is the offset strain tensor, and ejf and e°( are the axial and shear components 

respectively. Small offset strain values, on the order of 10 |i£, represent the initiation of 
yielding and are nearly the same as the proportional limit definition of yielding. 

However, larger offsets, such as the 0.2% offset criterion, give a macroscopic view of 
yielding and overall plastic flow. Because these large values cause considerable changes 
in material state the specimens cannot be used to determine multiple yield points in a 
locus. As a result, small offset strains are often used in yield surface experiments. Target 
offset values differ between experimental investigations. Some examples are. Helling et 
al. [26-28] used 5 p£, Gil et al. [29-30] used 30 pe, and Nouaihas and Cailletaud [31] 
used 100 |J£. 

The back extrapolation method is probably least used of the methods described 
because it requires a material with near linear hardening characteristics. In addition, it is 
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often necessary to strain the test specimen well into the plastic region in order to fit a 
straight line. As a result, the test specimens can not be used for multiple yield point 
determinations. This technique was used by Taylor and Quinney [32] to determine the 
multiaxial yielding behavior of copper, aluminum, and mild steel. 

The next question when developing a yield surface testing program is what 
control method should be used. The two most obvious choices are constant loading rate 
and constant strain rate. Studies performed by Phillips and Lu [25] and Wu and Yen [33] 
determined that there was no significant difference between stress-controlled and strain- 
controlled yield surfaces for pure aluminum specimens tested on a servohydraulic test 
machine. 

Another consideration for yield surface studies is the effect of strain rate on 
yielding. Ellis et al. [34] performed a study on dependence of probing rate on the small 
offset yielding behavior of type 316 stainless steel at room temperature. It was found that 
for strain rates between 100 and 500 |ie/min there was no significant change in the yield 
surface. These results agree with classical rate-independent plasticity. However, 
plasticity is always rate-dependant to some degree and the rate dependence typically 
increases at elevated temperatures. This is especially true at high stresses. Therefore, the 
strain rate used during probing will play an increased role at elevated temperature. 

1.4.2 The Yield Surface 

Investigations of the initial yield surfaces of metals date back to before the 
aforementioned work of Taylor and Quinney [32], Their research found that for copper 
and aluminum the von Mises yield criterion more accurately described initial yielding in 
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the axial-shear stress plane than the Tresca criterion. This conclusion has been confirmed 
by many other researches, such as, Phillips et al. [22], Liu [35], and Helling et al. [26]. 

The evolution of the yield surface, as a material is subjected to permanent 
deformation, is also of interest. What these subsequent yield surfaces look like depends 
on the definition of yielding used. If a large offset is used, such as the 2000 |!£ used by 
Hecker [36], the subsequent yield surface appears to be an isotropic expansion of the 
initial yield surface. However, if a small offset or proportional limit criterion is used the 
subsequent yield surface exhibits a combination of isotropic, kinematic, and distortional 
hardening. Experimental work done by Khan and Wang [37] provides a study of the 
effect of offset strains ranging from 200 to 2000 |ie on yield surface shape. 

In addition, yield surfaces can also exhibit substantial cross-effects. Positive 
cross-effects, an increase in the size of the yield surface normal to the prestrain direction, 
were found in studies performed by Phillips and Tang [23] and Wu and Yen [33]. In 
contrast, negative cross-effects (decrease in width normal to prestrain direction) were 
seen by Michno and Findley [21] in mild steels. A study by Williams and Svensson [38] 
showed that aluminum exhibited large cross-effects when subjected to torsion 
prestraining, but zero cross-effects after axial prestraining. 

When analyzing subsequent yield surfaces it is quickly realized that isotropic 
hardening, kinematic hardening, and cross effects cannot fully explain experimental 
results. Subsequent yield surfaces also show considerable amounts of distortional 
hardening. This is characterized be a decrease in curvature of the yield surface in the 
direction opposite to the direction of prestraining. The yield surface also continues to 
exhibit symmetry about the axis of prestraining. These conclusions have been found by 
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many researches including Phillips et al. [22], Phillips and Moon [24], Helling et al. [26], 
and Wu and Yeh [33]. 
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Shear Stress 


von Mises Yield Theory 

Tresca Yield Theory 



A xial Stress 

Figure 1.1: Relative difference between yield surfaces predicted by von Mises and 

Tresca yield criteria. 
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Shear Stress 



Initial Yield Surface 
Isotropic Hardening 
Kinematic Hardening 


Axial Stress 

Figure 1.2: The effects of isotropic and kinematic hardening on an initial yield surface. 
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Shear Stress 


Prestrain Direction 



A xial Stress 

Figure 1.4: Change in curvature of the backside of a yield surface due to an axial 

prestrain with added cross effects. 
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Shear Stress Shear Strain 



Axial Stress 

Figure 1.5: Difference between the direction of the strain rate and stress rate during a 
strain controlled nonproportional loading. 
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Chapter 2 

Experimental Methods 


The yield surface evolution experiments carried out for this study were complex 
and tedious. The following discussion concerns the methodology used in this study to 
perform yield surface evaluations and prestraining of the specimens. This chapter will 
discuss the materials, test specimen, test equipment, and test procedures used. The goal 
was to apply the same methodology to a wide range of materials in order to develop yield 
surface evolution equations. 

2.1 Materials 

This study consisted of testing three materials with significantly different 
compositions. The materials tested were Haynes 188, 316 stainless steel, and Inconel 
718. By testing such different materials it was hoped that differences in the way the yield 
surfaces evolve could be seen. 

2.1.1 Haynes 188 

Haynes 188 (HN 188) is a cobalt-based alloy, see Table 2.1 for the composition 
given by the manufacturer. Its primary uses include nozzle guide vanes, stator blades, 
and combustion liners in aero-propulsion systems. The material was tested in a solution 
annealed state. In this state the microstructure contains second phase carbides that 
provide an additional hardening component to the material response. The solution 
annealing process was conducted as follows, 


NAS A/CR— 2001-211162 


29 


Anneal at 1 175°C for one hour. 




• Cool to room temperature in Argon. 

The elastic moduli for HN 188 as found from the current study are shown in Table 2.2. 
All the values in Table 2.2 are averages over each specimen tested from the initial yield 
surface experiments for each material. 



Fe 

Ni 

Cr 

Co 

C 

S 

Mn 

Si 

P 

SS 316 

67.08 

10.2 

17.7 

0.19 

0.04 

0.013 

1.75 

0.57 

0.03 

HN 188 

1.17 

22.66 

22.11 

±38.7 

0.09 

0.002 

0.80 

0.35 

0.005 

IN 718 

±18.6 

53.58 

17.52 

0.39 

0.034 

0.002 

0.12 

0.07 

0.006 


Mo 

B 

Cu 

w 

Nb+Ta 

Ti 

Al 

La 

N 

SS 316 

2.08 

“ 

0.28 

- 

- 

- 

- 

- 

0.067 

HN 188 


0.003 

- 

14.06 

- 

- 

- 

0.052 


IN 718 

2.87 

0.004 

0.05 

- 

5.19 

0.95 

0.57 

- 

- 


Note: SS 316 heat A12210, HN 188 heat 188061714, IN 718 heat S618 


Table 2.1: Chemical Composition of Materials Tested. 



E (GPa) 

G (GPa) 

HN 188 

187.0 

71.6 

SS 316 

171.2 

62.9 

IN 718 

144.4 

54.8 


Table 2.2: Mechanical Properties of Materials Tested. (Elastic and Shear Modulus) 


2.1.2 316 Stainless Steel 

Type 316 stainless steel (SS 316) is a stable austenitic alloy with relatively low 
strength that is primarily used for pressure vessels and piping. The SS 316 specimens 
were again solution annealed, but with the following procedure. 
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• Anneal at 1065°C for 30 minutes. 

• Cool to 537°C in 9 minutes. 

• Furnace cool to room temperature. 

The composition of SS 316 is shown in Table 2.1 and the elastic moduli in Table 2.2. 

2.1.3 Inconel 718 

Inconel 71 8(IN 718) is a very high strength nickel-based alloy. It is mainly used 
in aeropropulsion applications for structural components such as disks, blades, and shafts. 
IN 7 18 was tested in an aged state. In order to reach this state the specimens were first 
solution annealed and subsequently aged. The thermal processing was conducted as 
follows. 

• Anneal at 1038°C for 1 hour. 

• Air cool to room temperature. 

• Age in vacuum at 720°C for 8 hours. 

• Cool to 620°C at 55°C/min and hold for 8 hours. 

• Cool in Argon to room temperature. 

This heat treatment produces y" precipitates that significantly increase the hardness of 
the material. Table 2.1 shows the chemical composition of IN 718 and Table 2.2 gives 
the elastic moduli. 

2.2 Test Specimens 

This section details the two different types of test specimen used. It also 
discusses the procedure for preparing a specimen to be tested. 
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2.2.1 Specimen Dimensions 

All of the test specimens used where of the thin-walled tube type. However, two 
different size specimens were used. The specimen sizes have no effect on the testing 
procedure other than the physical setup of the load frame grips and geometric constants 
used for interpreting the test data. Both HN 188 and SS 316 used the specimen type 
referred to as the 2-inch specimen due to its two-inch diameter grip section. IN 7 1 8 is 
the specimen referred to as the 1-inch specimen due to its one-inch diameter grip section. 
The exact dimensions of the 2-inch and 1-inch specimens are shown in Figure 2.1. 

2.2.2 Specimen Preparation 

Prior to testing each new specimen needed to be dimpled in order to accept the 
biaxial extensometer and fitted with thermocouples. Dimpling consisted of placing two 
small indentations in the gauge section of the specimen for the extensometer probes to 
rest. This was performed by using a dimpling device custom built at NASA GRC. The 
equipment consisted of a base that firmly held the test specimen in place with a screw 
used to lower the dimpling tool (supplied with the extensometer) to the specimen. 

Torque was then applied to the screw via a torque wrench in order to dimple the 
specimen. This setup allowed the same dimpling force to be applied to each specimen, 
precisely controlling the size of the dimples. Table 2.3 show the dimpling torques used 
for each material. 


Material 

Torque (N*m) 

HN 188 

2.03 

SS316 

1.58 

IN 718 

1.81 


Table 2.3: Dimpling torque. 
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The second step in preparing the test specimens was to spot weld thermocouples 
to the gauge section of the specimen. A total of three thermocouples were applied to 
each test specimen. The thermocouples were applied at the center and 12.5 mm above 
and below center on the gauge section. The center thermocouple was used to control the 
induction heating system. The remaining two thermocouples defined the temperature 
gradient across the gauge section. 

2.3 Test Equipment 

All experiments in this study were conducted on a computer controlled MTS 
biaxial servohydraulic test machine. The machine is pictured in Figure 2.2 and is located 
in the NASA GRC multiaxial fatigue lab. The maximum capacities of the load frame are 
±220,000 N axial loading and ±2,260 N*m twisting moment. The test specimen is held 
by water-cooled, hydraulically actuated grips. The top grip is fixed while the bottom grip 
is attached to a hydraulic actuator that is capable of independent vertical translation and 
rotation. Two MTS 458 analog controllers, one for axial motion and one for torsional 
motion, control the actuator. Further details about the biaxial test machine are given by 
Kalluri and Bonacuse [39]. 

Specimen heating was accomplished by using a closed-loop induction heating 
system as described by Ellis and Bartolotta [40]. The system consisted of an Ameritherm 
15-kW radio frequency induction heating unit and three adjustable, water-cooled copper 
coils surrounding the gauge section of the specimen. The copper coils could be raised or 
lowered independently along the gauge section of the test specimen in order to obtain an 
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acceptable temperature gradient (± 1 % of the absolute test temperature). See Figure 2.3 
for a photo of the copper coils around a l-inch specimen. 

In order to conduct yield surface studies, strain measurements accurate to the 
microstrain level are required. This is because loading must be quickly stopped once 
yield is detected in order to ensure the material state is not significantly disturbed. In 
addition, the device used to measure strain must also perform on a micro strain level at 
elevated temperatures. This depends not only on the strain measurement device, but 
also on the ability of the heating equipment to maintain a constant temperature and the 
elimination of electronic noise. 

This study utilized an MTS water-cooled biaxial extensometer capable of 
operating over a large temperature range. The extensometer used two alumina rods, 
spaced 25mm apart to measure axial deformation and twist. The alumina rods fit into 
the dimples placed in the test specimen as discussed earlier. Lissenden et al. [41] supply 
further details on the biaxial extensometer. The biaxial extensometer mounted to a 
specimen is shown in Figure 2.3. 

Custom written software and a PC were used to control all experiments. The PC 
was equipped with analog-to-digital (A/D) and digital-to-analog (D/A) hardware. The 
D/A hardware sent a command signal from the PC to the load frame while the A/D 
hardware received test data. Both sampled at 1000 Hz. The test data received by the 
A/D hardware averaged the data over every 100 points in order to minimize the effects 
of electronic noise, providing 10 data points per second written to output files. 
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The test program utilized two separate control programs. One was used to 
determined the yield points used to construct the yield loci while the other program 
preformed the preloadings. 

2.4 Strain Controlled Load Paths 

All three materials were subjected to a series of three strain controlled load paths. 
Path I is purely axial strain, Path II is purely shear strain, and Path HI is a non- 
proportional strain path. A schematic of each load path is shown in Figure 2.4. Along 
each load path several stops were made in order to conduct a pair of yield surface 
determinations. These points are indicated by letters in Figure 2.4. In addition, each load 
path was cycled either two or three times depending on the hardening characteristics of 
each material. For each load path the maximum equivalent strain is 15,000 (J£, shown by 
a dashed circle in Figure 2.4. 

Three additional load paths were carried out on IN 718. These paths were not 
carried out on the other materials due to time constraints and the lack of success with IN 
718. These additional load paths all subject the material to both axial and torsional strain. 
They are designated Path IV, Path V and Path VI and are shown in Figure 2.5. 

All six load paths used an equivalent strain rate of 100 |i£/s. For axial-torsional 
loading equivalent strain is given by. 



where £ IJ is the strain rate tensor, v is Poisson’s ratio, and £ u and £ ]2 are the tensorial 
axial and shear strain rates respectively. 


NASA/CR— 2001-211162 


35 


2.5 Yield Loci 


Each yield locus was constructed by probing for yielding in 16 unique directions 
in the axial-torsional stress plane. The angles used for probing were identical for each 
locus, however the order of probing varied depending on the prestrain history of the 
material. The angles were always setup in such a fashion that the first probe was in a 
direction normal to the prestrain direction. Furthermore, the order of subsequent probes 
was chosen in order to minimize the changes to the material state. For example. Figure 
2.5 shows the order of probing following a torsional prestrain. In this pattern, the angle 
between probes was either 180° or 90° in the hopes that the effects of the probes would 
counteract each other. Minimization of the changes in material state while performing 
yield surface determinations is further discussed by Hecker [42]. In addition, each yield 
surface was repeated to verify that the material state was undisturbed. 

Each of the 16 individual yield probes can be broken down into the following 
three step process, 

(i) Use a least squares fit to determine the coefficients of the elastic loading 
line over a predefined stress range within the elastic region. These 
coefficients include the elastic modulus, E, shear modulus, G, and initial 
stresses <7,°, and <7,° 2 . 

(ii) Continuously calculate the offset strain components using the following 
equations, 
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ff 12 - Cf [ f 


(2.3) 


-«/r 


= e i2 


2C 


From this determine the equivalent offset strain (Equation 1.33) and 
compare it to the target value of 30 |i£. 

(iii) Once the equivalent offset strain exceeds the target value unload the probe 
to the starting position. 

The 30 |i£ value for the offset strain was chosen based on previous experience 
from studies performed by Lissenden et al. [43] and Gil et al. [44]. In order to obtain 
optimum results with the least scatter each material used a different probing rate. Prior to 
testing a material a rate study was performed in order to determine the optimum probing 
rate. This was done by conducting initial yield surfaces tests for each material at various 
stress rates. Multiple runs of each rate were perform and the rate that showed the best 
repeatability was used as the yield surface probing rate for that material. The results of 
these rate studies are given in Table 2.4. 


Material 

Equivalent Stress Rate 
(MPa/s) 

Equivalent Elastic Strain 
Rate (um/m/s) 

IN 718 

2.07 

10 

HN 188 

17.9 

100 

SS 316 

7.24 

50 


Table 2.4: Yield surface probing stress rates. 
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Figure 2.2: Axial-Torsional test machine and MTS 458 analog controller. 
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Figure 2.3: Gripped 1-inch specimen with mounted biaxial extensometer and induction 

heating coils. 
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Figure 2.4: Strain controlled paths, (a) pure axial Path I, (b) pure shear Path 0, (c) 

combined axial-shear Path ID. 
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Probe 

Number 

Probe 

Angie 

Probe 

Number 

Probe 

Angie 

1 

12? 

9 

79° 

2 

192° 

10 

259" 

3 

102° 

11 

170“ 

4 

282° 

12 

350° 

5 

57° 

13 

125° 

8 

237° 

14 

305° 

7 

147° 

15 

35° 

8 

327° 

18 

215° 


Figure 2.5: Example of probing order after a torsional prestrain. 
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Chapter 3 
Data Analysis 


The overall goal of this research was to describe the evolution of the yield surface 
when subjected to plastic deformation. As a result, an integral part of the data analysis 
was to find a way to describe the yield surface by a model with parameters that change 
as the yield surface evolves. The first section of this chapter will discuss some of the 
attempts made to find such a model. Section 3.2 will discuss how the model developed 
by Voyiadjis et al. [13] and its parameters were used to reproduce yield surfaces and how 
the model parameters affect its shape, size, and position. Section 3.3 will detail the steps 
taken to reduce the model to an axial-torsional form that could be fit using linear 
regression. Finally, Section 3.4 will discuss the statistical methods used to quantify 
goodness of fit. 

3.1 Mathematical Representation of Yield Surfaces 

The first attempts made at fitting a shape to the yield surface data collected in 
this study applied the simplistic approach of fitting a circle with a flattened backside. 

The definition of the backside of a yield surface is the portion of the yield surface from 
approximately 90 to 270 degrees from the prestrain direction. This initial attempt to 
describe the yield surface was based on two principles; (i) the initial yield surface as well 
as the front side of the yield surface are circular in the Oi i-V3g 12 stress plane, (ii) the 
backside of the surface flattens compared to the front and can be represented by a straight 
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line. This representation proved adequate for describing yield surfaces with little or no 
distortion, however, it quickly faltered when exposed to highly distorted yield surfaces. 
In addition, the backside of the yield surface is not a straight line even when significant 
distortion is present. 

The second approach involved the regression of polynomials. Since the problem 
at hand is highly nonlinear the polynomials needed to be of order 2 or greater. 
Furthermore, terms involving both axial and shear stress needed to be included in order 
to obtain the distorted shape of the yield surface. As a result, the regression parameter 
became very complex. Nonetheless, the major shortcoming of the approach came from 
the fact that it was not always possible to obtain a continuous smooth surface. When 
used to predict a yield surface, the polynomial regression often produced a cusp when 
the curve intersected the axial-stress axis. This was not consistent with any of the 
experimental results collected. 

The polynomial approach was extended to include conformal mapping. It was 
hoped that a data conversion could be found that would allow the data to be fitted and 
then converted back to its original form. This process again produced non-continuous 
yield surfaces. In addition, the mapping process increased the error between data and 
fit curves resulting in poor yield surface representations. 

Finally it was decided that a yield function approach was needed. Several yield 
function based distortional models were discussed in Section 1.2.5. Of these, the model 
developed by Voyiadjis et al. [12, 13] was chosen for two reason: ( 1 ) it produced a 
continuous convex yield surface, (2) all of the pertinent yield surface characteristics 
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could can be controlled by four parameters. From this point this distortion model will be 
referred to as the Voyiadjis model. 

3.2 Yield Surface Representation 

The algorithm used to reproduce a yield surface once the model parameters were 
determined did not require a reduction to the axial-torsional form. The Voyiadjis model 
was used in the form shown in Equations 1.20 and 1.21 and programmed in to MathCAD. 
The yield surface was predicted by incrementing stress in the same probing direction 
used in the experimental portion of this study. First the program automatically found the 
starting stress point for all probes. To do this, probes 1 and 2 were fit with a line as well 
as probes 3 and 4. The intersection between these two lines gave the initial stress point 
for probing, as shown in Figure 3.1 where C represents the center of probing. From this 
point probing is carried out in all directions. This technique produced yield points in the 
same directions as the experimentally determined yield points. 

Each of the four parameters control a different type of hardening. The parameter 
a is related to isotropic hardening. Its value does not have a bound, but a-\ corresponds 
to a radius equal to the initial yield strength averaged from all specimens tested, k 0 . Then 
if a>\ the yield surface is smaller than the initial surface and if a< 1 the yield surface is 
larger than the initial surface. This is shown in Figure 3.2. The parameter d determines 
the amount of distortion there is in the backside of the yield surface. A value of zero 
gives no distortion. If the value is greater then zero than the yield surface distorts with an 
increase in curvature along the backside, as shown in Figure 3.3. If the value of d is less 
then zero the yield surface distorts with a decrease in curvature, as shown in Figure 3.4. 
The backstress parameters, an and a^, control the location of the center of the yield 
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surface. The values of these parameters depend on the properties of the material being 
tested, but they generally have an association with the maximum stress reached during 
a preloading. 

The value of the directional unit tensor, u,j, was calculated from the prestrain 
direction. Since all preloadings occurred in the axial-shear strain plane the axial and 
shear proportions of the loadings could be found from the cosine and sine of the loading 
angle respectively. If the directional tensor is given as, 

xcosd A-sinf? 0 

V 0 = a - sin 6 0 0 (3.1) 

0 0 0 

where 0 is the prestrain direction. By enforcing Huyll = 1 the value of x can be found 
from, 

A 2 cos 2 # + 2A 2 sin 2 # = l. (3.2) 

3.3 Parameter Determination for the Voyiadjis Model 

The first step in determining the parameters for the Voyiadjis model was to 
reduce it to its axial-torsional form by making the appropriate assumptions. Since the 
entire set of yield surface data is in the axial-torsional stress space the stress and 
backstress tensors can be reduced to, 

°u a i2 0 ] \a n a i2 0 

v,j = Or. 0 0 «ij = a, 2 0 0 . (3.3) 

0 0 oj [000 

The deviatoric stress and deviatoric backstress are. 
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One additional assumption is that the direction tensor, U;j, is also contained in the axial- 
torsional plane giving. 


v :: 


v, 


V t 2 

0 


v l2 

0 

0 


0 

0 

0 


( 3 . 5 ) 


When the above assumptions are made only 25 of the 81 components of the tensor My*/ 
(Equation 1 .21 ) are non- zero. Next, My*/ is substituted into Equation 1 .20 and the 
number of non-zero components furthers reduces to 1 1 for the completed axial-torsional 
form of the Voyiadjis model. 

With the axial-torsional reduction completed, work was initiated to obtain a form 
that could be used with a regression technique in order to obtain the model parameters. 
The problem was originally approached with the desire to obtain all four model 
parameters directly from regression. However, it is assumed that the backstress follows 
the same direction as the prestrain direction. As a result, it was easier to use an iterative 
method to find the backstress. An initial guess was input and the value was incremented 
in the direction of prestrain until the fit no longer improved. The method used to quantify 
the goodness of fit will be discussed in the next section. Since the backstress was not 
included in the regression the model reduced to a simple linear problem with the 
following characteristic equation, 
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Y = ft 0 k S +fi t Z 


( 3 . 6 ) 


where. 
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( 3 . 11 ) 


( 3 . 12 ) 


and ko is the initial yield strength averaged over all specimens of the same material tested. 
Linear regression was applied to Equation 3.6 to obtain following system of equations 
that can be solved by using matrix inversion, 


' n-k 2 XZljAoljin 

k 2 iz xz : JUj jxzzj 


( 3 . 12 ) 


where n denotes the number of probes in a yield surface. Finally Equations 3. 1 1 and 3.12 
are used to solve for a and d. 
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The yield surface fits were produced using both repeated yield surfaces at each 
point along a load path. This provided a maximum of 32 separate yield points in 16 
directions for each yield surface. Prior to, fitting the data were inspected visually to 
determine if a set contained outliers. A point was considered an outlier if it clearly did 
not fit the pattern suggested by the remaining data points. If a point was suspected of 
being an outlier fits were conducted with and without the data point and the better of the 
two was accepted. Yield points were excluded on an individual basis. If a yield point 
was excluded the corresponding point in the repeated surface was not necessarily 
excluded. 

3.4 Goodness of Fit Statistics 

The next problem was to quantify the goodness of the fit reproduced by the 
outlined process. Since linear regression was used the correlation coefficient was found 
to be sufficient for analyzing the fits. The experimental data and model predictions had 
to first be converted to fit the form given by Equation 3.6. This was accomplished by 
substituting the data and model parameters into Equations 3.7 and 3.12. The correlation 
coefficient statistic in terms of Y and Z is, 


R = 


£ZY- 


I Z£Y 


V(lZ 2 -nZ 2 X IT 2 -nE 2 ) 


( 3 . 12 ) 


where Z and Y are the average of all Z and Y values respectively. The value of R was 
then squared to get the coefficient of determination. It is stated by Kiemele et al. [45] 
that \R\ > 0.7 represents a good fit. However the yield surface data did not produce R' 
values that large, but still produced good fits in the axial-torsional stress plane. The 
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actual R 2 values range between 0.4 and 0.15. There is a correlation between R 2 and cycle 
number such that the R 2 values tend to decrease in the later cycles. This is largely due to 
an increase in experimental scatter in the later cycles. In addition, the fits used were not 
always statistically the best fit. This is due to the assumption that d< 0. All fits used were 
the best fit where this assumption was satisfied. 

3.5 Yield Surface Fitting Methodology 

Each yield surface was fit using the same series of steps. First an initial guess for 
the backstress parameters was made. This was accomplished by fitting a circle to the 
front side of the yield surface and using the center of the circle as the backstress. This 
worked very well for initial yield surfaces and for surfaces in the first cycle of a strain 
path. However, surfaces from later cycles sometimes required the initial backstress to be 
determined from experience. This was because as the distortion increased in later cycles 
the backstress no longer corresponded with the center of the front side. Therefore the 
initial backstress was obtain from the patterns observed in the evolution of the backstress 
from the initial cycle. 

Once an initial backstress was found the yield function was fit to the experimental 
data. The fitting process produced the isotropic parameter, a, and the distortional 
parameter, d , for a given backstress. Next, the backstress was incremented in the 
direction of prestraining. Because of this iterative process it was often desirable to make 
the initial guess for the backstress slightly under predict the magnitude of the expected 
final backstress value. This enabled the generation of data on both sides of the optimum 
fit. The best fit was the one with the largest R 2 statistic and a negative value of d. 
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Once the best fit was obtained, it was compared to the experimental data in order 
to search for possible outliers. An outlier was considered to be a yield point that did not 
follow the pattern given by the other data. Outliers were determined largely by judgment 
However, if a yield point was suspected it was removed from the fit and the entire 
process outlined above was conducted again. The new fit, without the outlier, was then 
compared to the original. The fit that had the largest R~ statistic was then chosen as the 
fit for that yield surface. 


N AS A7CR— 2001-21 1 162 


53 



# 



- 400.00 0.00 400.00 

Axial Stress (MPa) 


Figure 3.1s HN 188 Path n, Point C, Cycle 1 showing how the intersection of two lines 
between yield points gives the center of probing. 
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Figure 3.3: Distortion of a yield surface with changing d ranging from 0.0 to -0.2. 

(ko = 203. 1 MPa) 
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Figure 3.4: Distortion of a yield surface with changing d ranging from 0.0 to 0.01 

(k 0 = 84.4 MPa) 
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Chapter 4 

Results and Discussion 


With the amount of experimental data collected it is important that the results be 
presented in a concise manner. It was decided that the data should be given in a format 
that handled each strain path separately. This allows for comparison between the three 
materials for each strain path. Section 4. 1 will first introduce some of the experimental 
yield surface data. Each section from there will briefly discuss the behavior of each 
material for a given load path and then show summary plots of the evolution of the 
primary variables. The variables of interest are a, d, an, and an because these are the 
four parameters that describe the yield surface in the Voyiadjis distortion model. After 
the three primary strain paths are discussed the additional paths applied to IN 7 18 will be 
discussed. Note, this section does not try to present all of the yield surface data collected. 
An entire set of the data and model fits can be found in Appendix A. 

4.1 Yield Surface Data 

All yield surface tests revolve around the ability to obtain individual yield points 
from each probing direction. Therefore, it is appropriate to show details of at least one 
yield probe. Figure 4.1 shows (a) the stress-strain response and (b), the offset strain-total 
strain response of a 12° yield probe for HN 188. From Figure 4.1 (a) the small offset 
strain is visible from the separation between the loading and unloading lines of the stress- 
strain curve. Furthermore, Figure 4. 1 (b) shows that the offset strain reached a value of 
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31.6, slightly overshooting the target of 30. This discrepancy is acceptable and does not 
constitute a large deviation from the target. In addition. Figure 4.1 (b) shows that the 
offset value slightly decreases during the unloading of the probe. This decrease can be 
contributed to vicoplasticity effects due to the elevated temperature of the testing. It 
could also be the result of a slightly nonlinear elastic stress-strain response. 

The first step for all yield surface evolution tests was determining the initial yield 
surface. The initial yield surface was assumed to be circular in the o ]r V3*Oi2 stress 
plane. Therefore, each initial yield surface was regressed to a circle by setting the 
directional unit tensor, in the Voyiadjis model equal to zero. This caused the 
distortion portion of the model to go to zero and produce a circular yield surface. Figure 
4.2 shows an initial yield surface and model prediction for SS316, Path I. In this figure 
there is one yield point that falls well inside of the fitted surface. This is typical of the 
experimental scatter inherent to this study. 

Each of the three materials has considerably different initial yield surfaces. All 
three have different radii and the initial yield surface for IN 718 is translated in the 
compressive direction due to its strength differential. More information on the strength 
differential in IN 7 1 8 and other materials can be found in Gil [46] and Hirth and Cohen 
[46]. Figure 4.3 presents the average yield surface fit for all of the specimens tested. 

The difference in initial yield strengths and the strength differential are clearly visible. 

Different types of subsequent yield surfaces were produced for each of the six 
prestrain directions included in paths 1-ID. Figures 4.4-4.6 shows subsequent yield 
surfaces for each of these prestrain directions for HN 188, SS 316, and IN 718, 
respectively. 
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The solid line in each figure represents the initial yield surface. From the graphs 
it is clear that the kinematic hardening component of the yield surfaces follow the 
prestrain direction. It can also be seen that the distortion is in the opposite direction of 
the prestrain. In addition, when compared to the initial surface each subsequent surface 
shows isotropic hardening. Closer inspection of the data shows that each yield surface 
has two distinct parts: ( 1 ) a well-defined circular front side, and (2) a poorly defined 
flattened backside. The front side of the surfaces encompasses the yield point within the 
range of ±90° from the direction of the most recent prestrain direction. The front side is 
also more-or-less centered with respect to the prestrain direction. Along the backside 
of the yield surface (approximately 90° to 270°) the curvature decreases. Significant 
experimental scatter in the yield points along the backside of the yield surface was 
typically observed. This scatter tended to increase as the specimen accumulated plastic 
strain as it proceeded through a given load path. 

4.2 Effect of Load Path Cycling 

Each load path was cycled for two or three cycles. The number of cycles was 
determined by when the yield surfaces stopped evolving. Each material behaved 
differently when the load paths were cycled, but in general, HN 1 88 and IN 7 1 8 required 
2 cycles and the SS 316 required three cycles for cyclic hardening to stop. In fact, IN 718 
stopped evolving by the end of the first cycle, but a second was conducted to confirm 
this. Figure 4.7 shows the cyclic evolution of the yield surfaces at selected points for HN 
188, Figure 4.8 for SS 316, and Figure 4.9 for IN 718. From these figures it can be seen 
that the yield surfaces changed the most in the first cycle. Each consecutive cycle then 
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produced less hardening effects. SS 316 continued to harden through all three cycles. 
HN 188 showed considerable hardening between the first and second cycles, but little 
between the second and third. As a result, only Path II was carried out for three cycles 
for HY 188. Path II and Path HI were stopped after just two cycles. The hardening of 
IN718 nearly stopped after just one cycle. The yield surfaces for cycles two virtually 
overlap those for cycle one. In general, the number of cycles needed for the yield surface 
evolution to stop was inversely related to the initial yield strength of the material. 

4.3 Strain Path I 

Strain path I was a purely axial strain path as shown in Figure 2.4 (a). Both HN 
188 and IN 718 were subjected to 2 cycles of this path while SS 316 was cycled three 
times. In addition, time constraints made it necessary to conduct the yield surface 
determination for HN 1 88 and IN 7 1 8 at only points A, C, E, G, and I. (see Figure 2.4a). 
Figure 4. 10 (a) shows the axial strain paths produced by Path I. Only one example for 
each load path is shown in Figure 4. 10 because the strain paths were well controlled and 
the results for all materials look identical. All cycles are shown and overlap. Figure 4.1 1 
provides the axial stress-axial strain response of each material for Path I. It clearly shows 
the relative amount of isotropic hardening each material exhibits. Both HY 188 and SS 
316 show a continuous increase in the maximum stress while IN 718 quickly reaches a 
maximum value. Of primary concern in this section is the evolution of the parameters 
for the Voyiadjis model. The parameters are shown in Figure 4.12-4.14. The data is 
presented in the order of the yield surface determinations along the load path. From 
Figures 4. 13 it can be seen that SS 316 has the most isotropic hardening because it has 
the largest change in the parameter a. SS 316 is also the only material for which a 
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becomes negative. HN 188 has the second largest change in a and IN 718 the least. For 
each of the three materials a decreases the most during the first cycle and then begins to 
level off. This mirrors the observations made about the cyclic evolution of the materials 
in section 4.2. 

The parameter d does not have a well-defined trend. Since this term measures the 
distortion in the yield surface it is greatly influenced by the scatter in the data for the 
backsides of the yield surfaces. The most well behaved test, from the standpoint of the 
distortion parameter, d , was HY 188 Path 1 (Figure 4.12). It shows d as nearly constantly 
decreasing. This indicates an increase in distortion along the backside as the load path 
progressed. The other two materials show this through roughly the first cycle but then it 
becomes difficult to see a trend. Most tests included a few yield surfaces that had more 
distortion than those around them. These yield surfaces often showed considerable 
experimental scatter and outliers that were not used in the fit. 

The axial backstress component, an, follows the same pattern as the axial stress. 

It starts at a value of zero for the initial yield surface for both HY 188 and SS 316. 
However, due to the strength differential shown by IN 719 an has a value of -145 MPa 
for the initial yield surface. After the initial yield surface, an takes on the same sign as 
the maximum axial stress obtained during the prestrain prior to a given subsequent yield 
surface. In addition, ai i increased in magnitude with each consecutive prestrain in the 
same direction. The torsional backstress component, a^, is not shown for Path I because 
it does not evolve during the load path. 

It is worth discussing here that IN 718 showed considerably more data scatter 
than the other materials. The additional scatter did affect the quality of the fits that were 
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obtained. For some of the IN 718 yield surfaces the front side contained scatter as well 
as the backside. Several attempts were made during testing to eliminate this problem. 
The yield surfaces for IN 7 1 8 were run at considerably lower stress rates than the other 
materials (see Table 2.4 for rates). This did reduce the scatter somewhat, but it is present 
in all IN 7 18 yield surfaces. 

4.4 Strain Path II 

Strain Path II is pure shear and is shown in Figure 2.4 (b). IN 718 was subjected 
to two cycles of this strain path while HN 188 and SS 316 were tested for three. Again, 
IN 7 18 was tested only at every other point shown in Figure 2.4 (b). Figure 4. 10 (b) 
shows an example of the axial strain-shear strain response for this load path. The shear 
stress-shear strain response for each material is shown in Figure 4.15. Figures 4.16-4.18 
show the evolution graphs for the Voyiadjis model parameters for HY 188, SS 316, and 
IN 7 18 respectively. 

In general, the same observations for the evolution of the parameters a and d were 
found for Path II as Path I. Parameter a again decreased as the strain path proceeded with 
the same relative changes between materials as seen with Path I. The evolution of a for 
HN 1 88 showed more scatter than for Path 1 . The distortion parameter again provided no 
definite trend. It again decreased at the start of the load path and then became scattered 
in later cycles. 

Path II involved the evolution of the torsional backstress component, a 12 . It 
followed the shear stress as ai i followed the axial stress for Path I. The axial backstress 
component did not evolve during this load path, but for IN 718 it had a non-zero value for 
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the entire test. This is again a manifestation of the strength differential of IN 7 18. The 
axial backstress maintained a constant value of -145 MPa. 

4.5 Strain Path III 

Unlike Path I and Path n. Path in subjects the test specimen to a nonproportional 
loading. Points B and F follow prestrains at 45 degrees and point D follows a prestrain at 
315 degrees (see Figure 2.4(c)). The axial strain-shear strain response is shown in Figure 
4.10 (c). Both the axial stress-shear strain and shear stress-shear strain graphs are shown 
for each material in Figure 4. 1 9. Figures 4.20-4.24 show the evolution graphs for the 
Voyiadjis model parameters for HY 188, SS 316, and IN 718 respectively. Note that 
Figures 4.20 and 4.21 show repeats runs of Path HI for HN 188 on separate specimens. 

In addition, Figure 4.22 shown the two runs of HN 188 overlaid on the same graph. 

All four model parameters evolved during this load path. Both a and d showed 
the same response as with the other load paths. The model parameter a decreased while d 
exhibited scatter and showed no real trend. The scatter in d was particularly evident with 
the yield surfaces at points C and E. The distortion parameter grew large at these two 
points. The yield surfaces for Path HI showed more experimental scatter than the other 
load paths. This was due to the fact that Path HI imparted more plastic strain on the test 
specimen per cycle than the other load paths. 

The backstress terms both evolved for this load path. Again they closely followed 
the movement of their stress counterparts. The axial backstress followed the axial stress 
and the torsional backstress followed the torsional stress, both at a lower magnitude but 
with the same sign. 
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Unfortunately the test of Path III on IN 7 18 was not completed for even one cycle. 
The test specimen buckled during the prestrain between point D and E. This is likely due 
to a difference in the heat treatment of the specimen used for the test. The specimen was 
aged separately from the test specimens used for the other IN 718 tests. The same 
procedure was followed, but at a different time. For an unknown reason the specimen 
had a yield strength, in the virgin state, that was higher than the other IN 718 specimens. 
It is suspected that this increase in initial yield strength also decreased the ductility of the 
specimen. 

4.6 Additional Inconel 718 Strain Paths 

For IN 7 18 three additional load paths were attempted. They were created as an 
attempt to gain additional data from specimens needed for other studies. See Figure 2.5 
for schematics of these additional strain paths. Only Path IV was successfully carried out 
for an entire cycle. The other two strain paths caused the specimens to buckle during 
prestraining. The evolution of the Voyiadjis model parameters for Path IV is shown in 
Figure 4.25. 

4.7 Summary of Experimental Results 

The materials tests provide a wide range of microstructures, compositions, and 
material properties. However, the evolution of the yield surface of all three materials can 
be described by a combination of isotropic, kinematic, and distortional hardening. 
Therefore, it appears that the same type of hardening laws can be used to describe each 
alloy. This study also has shown that distortional hardening is difficult to accurately 
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capture experimentally. The scatter in the backsides of the yield surfaces was relatively 
small for the first cycle, but became significant in later cycles. This makes it difficult to 
discuss the evolution of distortional hardening parameter with any confidence. The 
remaining parameters, a, an, and a.\ 2 , are all well behaved and there is high level of 
confidence in the results. 

4.8 Parameter Evolution 

Of the four yield surface parameters in the Voyiadjis model three have regular 
patterns to their evolution. These are a, a\ 1 , and a[ 2 . As a result, it was attempted to find 
equations to represent this evolution. 

This first parameter that was considered was a. It is clear that a decreases as the 
material is subjected to more plastic deformation. Therefore, it is logical to assume that it 
has a relationship to the plastic work where plastic work is defined as, 

W =<j£ p . (4.1) 

p ij lj 

When a was plotted versus plastic work the relationship proved to be of an exponential 
nature. Because a is not always positive there was a problem fitting an exponential 
equation to the data. To eliminate the problem the data was shifted by adding the 
absolute value of the largest negative value of a to the data. Figures 4.26-4.28 show the 
results of exponential fits for each material. For each fit all of the experimental data from 
each load path was used. The equations of the lines and R 2 statistics for each material are 
as follows, 

1. HN 188 

a = 1 . 1 85 * e-°- 04tHWp R 2 = 0.42 1 (4.2) 
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0 


SS 316 


a = 1 .021 * e^ omWp R 2 = 0.465 (4.3) 

3. IN 718 

a = 0.7664 * e ~° m ™ Vp R 2 = 0.7340 (4.4) 


The back stress parameters were more difficult to fit to evolution equations. This 
was because during each load path the sign of the backstress alternated several times. As 
a result it was necessary to use the sign function to define the evolution of the backstress. 
The sign of the stress at a given point was taken while the magnitude of the backstress 
was related to the equivalent total strain. This relationship was used because the sign of 
the backstress was always the same as the stress, but its magnitude increased as total 
strain increased. Figures 4.29 and 4.30 show the results of this procedure for HN 1 88 and 
the evolution equations obtained are, 

a„ = 1 13.18 + 3.19Z-0.035Z 2 +0.00084Z 3 + 3.2210‘ 6 Z 4 (4.5) 

a 12 = -67.53 + 3.33F + 0.89F 2 -0.001 IF 2 -7.86 10~ 6 F 4 (4.6) 

where. 


Z — (7 ll £ n 
F ^ [2 

The R 2 statistics for Equation 4.5 and 4.6 are 0.79 and 0.33 respectively. 


(4.7) 
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Figure 4.1: Stress-Strain, (a), and Offset Strain-Axial Strain, (b), response for an 

individual yield probe. 
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Figure 4.2: Initial yield surface for SS 316 showing model fit and parameters. 
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Figure 4.3: Initial yields surfaces for HN 188, SS 316, IN 718. 
(Averaged over all specimens tested.) 
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Figure 4.4: Examples of HN 188 yield surfaces. 
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Figure 4.7: Cyclic evolution of HN 188 at selected points. 
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Figure 4.8: Cyclic evolution of SS 316 at selected points. 
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Figure 4.9: Cyclic evolution of IN 7 1 8 at selected points 
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e 4 . 11 : Axial Stress-Axial Strain response for Path I. 
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Figure 4.12: Voyiadjis model parameter evolution for HN 188, Path I, HYII 89. 
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Figure 4.13: Voyiadjis model parameter evolution for SS 316, Path I, 610-01. 
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Figure 4.14: Voyiadjis model parameter evolution for IN 718, Path I, IN- 16. 
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Figure 4.16: Voyiadjis model parameter evolution for HN 188, Path II, HYE-86. 
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Figure 4.17: Voyiadjis model parameter evolution for SS 316, Path II, 610-04. 
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Figure 4.18: Voyiadjis model parameter evolution for IN 718, Path n, IN-23. 
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Figure 4.20: Voyiadjis model parameter evolution for HN 188, Path m, HYII-90. 
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Figure 4.21: Voyiadjis model parameter evolution for HN 188, Path HI, HYII-82. 
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Figure 4,22: Voyiadjis model parameter evolution for HN 188, Path HI for runs 1 and 2. 
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Figure 4.23: Voyiadjis model parameter evolution for SS 316, Path HI, 610-05. 
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Figure 4.24: Voyiadjis model parameter evolution for IN 718, Path HI, IN-27. 
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Figure 4.25: Voyiadjis model parameter evolution for IN 718, Path IV, IN-11. 
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Figure 4.26: Evolution of a versus plastic work for HN 188. 


NASA/CR— 2001-211162 


94 




1.00 



0.00 20.00 40.00 60.00 80.00 100.00 

Accumulated Plastic Work 

Figure 4.28: Evolution of a versus plastic work for IN 718. 
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Figure 4.29: Evolution ofay/ versus stress times total strain squared for HN 188 
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Figure 4.30: Evolution of a l2 versus stress times total strain squared for HN 188 
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Chapter 5 


Conclusions and Future Work 


This chapter provides a summary of the major conclusions from both the 
experimental and analysis portions of this study as well as some suggestions for future 
work. 

5.1 Conclusions 

This study had two very distinct parts. First was the experimental portion and 
the other was the data analysis portion. Both has it own set of conclusions. 

Experimental 

• Plastic flow in Haynes 188, 316 Stainless Steel, and Inconel 718 
complicated and highly alloy dependent, with each material exhibiting 
varying degrees of each of the three components of hardening; isotropic, 
kinematic, and distortional. 

• All subsequent yield surfaces had two distinct parts; a well defined front 
side and a poorly defined backside. The front side was nearly circular in 
the axial-shear stress plane and the backside was flattened. 

• Subsequent yield surfaces continued to evolve through multiple cycles of 
a given load path. 
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• Even thought the three alloys tested were vastly different it appears that 
the same type of evolution equations could represent the multiaxial 
hardening of each material. 

Data Analysis 

• The Voyiadjis model adequately describes subsequent yield surfaces with 
only four parameters; one isotropic, one distortional, and two backstress 
parameters. 

• The isotropic parameter, a, and kinematic parameters, a\\ and a 12, evolve 
smoothly though the entire set of yield surface tests. 

• The distortional hardening parameter, d, is difficult to describe, partly 
because of the large amount of scatter in the yield points on the backside 
of the yield surface. 

• Evolution of the isotropic hardening parameter and kinematic hardening 
parameters can be reasonably well described by simple equations. 

5.2 Future Work 

The work presented here attempts to quantify the evolution of subsequent yield 
surfaces. This process includes many experimental and mathematically difficulties. 
However, now that the groundwork is set for determining the necessary parameters for 
quantifying yield surfaces there are many areas that require more attention. Of primary 
concern is improving the yield surface fits. Some suggestions are, 

• Reducing the amount of experimental scatter in the backside of subsequent yield 
surfaces. Two approaches would apply, 

1 . Study the effect of probing rate on backside scatter. 
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2. Study the effect of the offset strain target value on backside scatter. 

In either case the sensitivity studies would be carried out for subsequent yield 
surfaces instead of initial yield surfaces as done with the current study. 

• Further investigate the increased yield surface scatter observed in IN 7 1 8 as 
instead of SS 16 and HY 188. 

In addition to improving the experimental data further work is needed in 
developing evolution equations for the yield surface parameters. The ultimate goal is to 
develop evolution equations to be used in a viscoplasticity model to predict yield surfaces 
for a given plastic strain value. 
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Appendix A 

Catalog of Yield Surface Data 


Appendix A begins with a list of all specimens tested and the dimensions of each. 
The second section contains all yields surfaces collected along with the fit data for each. 

A.l Specimen Dimensions 


Material 

Specimen 

Number 

Test 

Path 

Outside 
Diameter (in) 

Inside 

Diameter (in) 

Haynes 1 88 

HYII 82 

III 

1.0241 

0.8669 


HYII 85 

1 

1 .0227 

0.8660 

HYII 86 

II 

1 .0238 

0.8656 

HYII 89 

1 

1 .0231 

0.8663 

HYII 90 

III 

1.0237 

0.8653 





316 Stainless 

610-01 

1 

1.0257 

0.8675 


610-04 

II 

1.0252 

0.8678 

610-05 

III 

1.0257 

0.8680 





Inconel 718 

11 

IV 

0.8225 

0.6244 


14 

VI 

0.8220 

0.6238 

16 

1 

0.8195 

0.6245 

23 

II 

0.8230 

0.6244 

27 

III 

0.8190 

0.6236 

28 

V 

0.8203 

0.6242 


Table A.l: Specimen dimensions and test matrix. 
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A.2 Yield Surface Data 


All yield surface data are presented here. Figures A.l through A. 12 show HN 
188 data. Figures A. 13 through A.27 show SS 316 data, and Figures A.28 through A.32 
show IN 718 data. The experimental data is represented by plus signs and diamonds 
(the first and second repeated surfaces), the fits are shown by a solid line, and the initial 
surface is shown by a dashed line. In addition, the model parameters a and d as well as 
the R‘ values are shown on each graph. 
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Figure A.2: Haynes 188, Path I, Point C-I, Cycle 2, HYEI-89. 
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Figure A.3: Haynes 188, Path El, Points A-E, Cycle 1, HYII-86 
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Figure A.4: Haynes 188, Path n. Points F-H, Cycle 1, HYII-86. 
(The data for point I was bad therefore not fit was performed.) 
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Figure A.5: Haynes 188, Path H, Points B-E, Cycle 2, HYH-86 
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Figure A.6: Haynes 188, Path n. Points F-I, Cycle 2, HYE-86 
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Figure A.7: Haynes 188, Path II, Points B-E, Cycle 3, HYII-86 
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Figure A.8: Haynes 188, Path n. Points F-I, Cycle 3, HYII-86 
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Figure A.ll: Haynes 188, Path HI, Points A-F, Cycle 1, HYTI-82 
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Figure A. 13: 316 Stainless Steel, Path I, Points A-E, Cycle 1, 610-01 . 
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Figure A.14: 316 Stainless Steel, Path I, Points F-I, Cycle 1, 610-01. 
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Figure A. 16: 316 Stainless Steel, Path I, Points F-I, Cycle 2, 610-01. 


NAS A/CR— 2001-211162 


126 



Shear Stress * I 732 ( M Pa ► Shear Stress * 1 .732 (MPa) 



-400-00 -200 00 0.00 200.00 400.00 -40000 - 20000 0 00 200.00 

Axial Stress <MPa> Axial Stress i M Pa) 



-400.00 -200.00 0.00 200.00 400 00 -400,00 -200.00 0 00 200.00 

Axial Stress (MPai Axial Stress (M Pa) 


Figure A.17: 316 Stainless Steel, Path I, Points B-E, Cycle 3, 610-01. 
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Figure A.18: 316 Stainless Steel, Path I, Points F-I, Cycle 3, 610-01 
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Figure A.19: 316 Stainless Steel, Path n, Points A-E, Cycle 1, 610-04 
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Figure A. 20 : 316 Stainless Steel, Path n. Points F-I, Cycle 1 , 610-04 
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Figure A.21: 316 Stainless Steel, Path n. Points B-E, Cycle 2, 610-04 
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Figure A.22: 316 Stainless Steel, Path n. Points F-I, Cycle 2, 610-04 
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Figure A.23: 316 Stainless Steel, Path n. Points B-E, Cycle 3, 610-04 
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Figure A.24: 316 Stainless Steel, Path H, Points F-I, Cycle 3, 610-04. 
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Figure A.25: 316 Stainless Steel, Path m, Points A-F, Cycle 1, 610-05 
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Figure A.26: 316 Stainless Steel, Path IQ, Points B-F, Cycle 2, 610-05. 
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Figure A.27: 316 Stainless Steel, Path HI, Points B-F, Cycle 3, 610-05. 
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Figure A30: Inconel 718, Path II, Points A-I, Cycle 1, IN-23. 
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Figure A.32: Inconel 718, Path IE, Points A-D, Cycle 1. 
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Appendix B 


Yield Surface Fit Data 


Cycle 

Point 

<y (MPa) 

e (m/m) 

a 

d 

all 

a12 

1 

A 

0 

0 

0.8077 

0 

0 

0 

1 

C 

314 

0,015 

0.6428 

-6.30E-05 

49 

-8 

1 

E 

-348 

-0.0013 

0.4289 

-1 .60E-04 

-94 

-10 

1 

G 

-422 

-0.015 

0.2691 

-4.10E-05 

-120 

-8 

1 

1 

406 

0.0013 

0.2308 

-8.90E-05 

73 

-14 

2 

C 



0.202 

mm 

101 

-7 

2 

E 


-0.0013 

0.1663 

-1.00E-04 

-131 

-5 

2 

G 

-489 

-0.015 

0.1238 

-9.80E-05 

-176 

-11 

2 

1 

458 

0.0012 

0.058 

-3.21 E-04 

90 

-10 


Table B.l: Haynes 188, Path I, HYD-89. 


- 2001-211162 
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Table B.2: Haynes 188, Path n, HYII-86. 


NAS A/CR— 2001-211162 


146 


















































































































































Cycle 

Point 

o (MPa) 

e (m/m) 

x (MPa) 

Y (m/m) 

a 

d 

all 

a12 

1 

A 

0 

0 

0 

0 

1.086 

0 

0 

0 

1 

B 

225 

0.0106 

105 

0.0077 

0.7768 

-3.90E-04 


14.7 

1 

C 


-0.0106 

37 

0.0077 

0.4822 

-3.63E-03 

-78 

2 

1 

D 





0.3152 

-1.16E-04 


-56 

1 

E 


-0.0106 

-32 

-0.0076 

0.3105 

-3.60E-03 


-7 

1 

F 


0.0008 

202 

0.0006 





2 

B 

358 

0.0106 

196 

0.0077 

0.1219 

-4.55E-04 

56 

48 

2 

C 

-507 

-0.0106 

35 

0.0077 

0.1146 

-6.27E-04 

-170 

0 

2 

D 

389 

0.0106 

-240 

-0.0076 

0.228 

-9.46E-05 

112 

-115 

2 

E 

-572 

-0.0106 

-40 

-0.0076 

0.098 

-4.87E-03 

-210 

-20 

2 

F 

376 

0.0008 

254 

0.0006 

-0.076 

-1.16E-03 

70 

75 


Table B.3: Haynes 188, Path El (First Run), HYE-90. 


Cycle 

Point 

a (MPa) 

e (m/m) 

X (MPa) 

Y (m/m) 

a 

d 

all 

a12 

1 

A 

0 

0 

0 

0 

0.962 

0 

-10 

0 

1 

B 


0.0106 

106 

0.0077 

0.7558 

-1 .25E-03 


12 

1 

C 


-0.0106 

21 

0.0077 

0.4358 

-1 .06E-03 

-72 

0 

1 

D 

314 

0.0106 

-168 

-0.0076 

0.2746 

-1.09E-03 

62 

-59 

1 

E 

-446 


-23 

-0.0076 

0.1785 

-2.66E-03 

-111 

-10 

1 

F 



209 

0.0006 

0.081 

-4.00E-04 

86 

73 

2 

B 


0.0106 

208 

0.0077 

0.053 

-5.38E-04 

97 

77 

2 

C 


-0.0106 

33 

0.0077 

0.0283 

-1.15E-03 

-168 

10 


D 

416 

0.0106 

-246 

-0.0076 

-0.038 

-9.86E-04 

100 


2 

E 

-589 

-0.0106 

-45 

-0.0076 

0.081 

-3.74E-03 


-60 

2 

F 

383 

0.0008 

270 

0.0006 

-0.109 

-1.37E-03 

84 

87 


Table B.4: Haynes 188, Path IE (Second Run), HYE-82. 
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Cycle 


o (MPa) 

e (m/m) 

a 

d 

all 

a12 

1 

A 


0 

0.8359 

0 

-4 

0 

1 

B 


0.0075 

0.7042 

-3.49E-04 


0 

1 

C 

150 

0.015 

0.482 

-9.60E-04 

29 

1 

1 

D 


0.0075 

0.3321 

-6.01 E-04 

-48 

-6 

1 

E 

-195 

-0.0006 

0.2986 

-1.03E-03 

-58 

0 

1 

F 

-209 




-64 

-4 

1 

G 

-221 

-0.015 

0.083 

-1.40E-03 


-8 

1 

H 

206 



-1.60E-03 

40 

-8 

1 

1 

216 




65 

-2 

2 

B 

231 

0.0075 

0.094 

-1.49E-03 

71 

-2 

2 

C 

244 


0.093 

-2.96E-03 

90 

-1 

2 

D 


0.0075 

0.004 

-2.00E-03 


-4 

2 

E 


-0.0006 

0.1101 

-2.70E-03 


-ii 

2 

F 


-0.0075 

-0.0026 

-1.69E-03 


4 

2 

G 


-0.015 

-0.036 

-1.40E-02 

-120 

-11 

2 

H 

245 

-0.0075 

-0.114 

-1.00E-03 

40 

0 

2 

1 

269 

0.0006 

-0.1 

-2.00E-03 

95 

-4 

3 

B 

275 

0.0075 

-0.0905 

-1.47E-03 

100 

-4 

3 

C 

283 

0.015 

-0.122 

-1.02E-02 


-2 

3 

D 

-260 

0.0075 

-0.123 

-8.50E-04 

-83 

-2 

3 

E 


-0.0006 

-0.141 

-2.37E-03 

-106 

-3 

3 

F 


-0.0075 

-0.131 

-2.69E-03 

-120 

1 

3 

G 

-308 

-0.015 

-0.15 

-8.05E-03 

-125 

0 

3 

H 

266 

-0.0075 

-0.185 

-4.31 E-03 

64 

0 

3 

1 

296 

0.0006 

-0.13 

-3.77E-03 

114 

0 


Table B.5: 316 Stainless Steel, Path I, 610-01. 
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Cycle 

Point 

IBB 

IBB 

II 

d 

all 

msm 

1 

A 

0 

0 

IBB 


0 

■OH 

1 

B 

70 

5.40E-03 

II 

-3.28E-03 

1 

ii 

1 

C 

75 

1.08E-02 

0.8603 

-5.02E-03 

0 

17 

1 

D 

-77 

5.40E-03 

0.6438 

-5.35E-03 

-7 

-15 

1 

E 

-87 

-4.00E-03 

0.5983 

-6.82E-03 

-1 

-23 

1 

F 

-96 

-5.40E-03 

0.4182 

-6.24E-03 

0 

-26 

1 

G 

-100 

-1.08E-02 

0.3328 

-6.16E-03 

-3 

-28 

1 

H 

98 

-5.80E-03 

0.188 

-2.84E-03 

5 

14 

1 

1 

101 

3.00E-04 

0.206 

-3.15E-03 

5 

26 

2 

B 

110 

5.40E-03 

0.131 

-2.76E-03 

0 

28 

2 

C 

120 

1.08E-02 

0.1743 

-2.59E-03 

0 

34 

2 

D 

-117 

5.40E-03 


-2.45E-03 

-2 

-31 

2 

E 

-124 

^Tooe-03 

0.844 

-3.1 IE-03 

0 

-35 

2 

F 

-126 

-5.40E-03 

0.0587 

-2.71 E-03 

0 

IPPPM 

2 

G 

-135 

-1 .08E-02 

0.151 

-3.31 E-03 

0 

H 

2 

H 

129 

-5.80E-03 

0.0087 

-2.75E-03 

0 

30 

2 

1 

131 




-1 

43 

3 

B 

138 

5.40E-03 

-0.046 

-2.80E-03 

2 

47 

3 

C 

150 


E&sai 


0 

51 

3 

D 

-140 

5.40E-03 

Klffll 


-10 

-42 

3 

E 

-150 

-4.00E-03 

mw$m\ 


-10 

-39 

3 

F 

-157 




-7 

-57 

3 

G 

-leiii 


-0.064 

-2.73E-03 

-7 

-62 

3 

H 

149 

-5.80E-03 

-0.163 

-2.44E-03 

-14 

32 

3 

1 

164 

3.00E-04 

-0.126 

-2.69E-03 

-12 

53 


Table B.6: 316 Stainless Steel, Path II, 610-04. 
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Cycle 

Point 

a (MPa) 

e (m/m) 

x (MPa) 

y (m/m) 

a 

d 

all 

a12 

1 

A 

0 

0 

0 

0 

0.9918 

0 

-5 

2 

1 

B 

101 

0.0106 

46 




22 

11 

1 

C 

-175 

-0.0106 

-9 

0.0076 

0.2841 

-4.31 E-02 

-51 

-5 

1 

D 

175 

0.0106 

-85 

-0.0076 

-0.0027 

-1.73E-04 

35 

-34 

1 

E 

-271 

-0.0106 

-3 

-0.0076 

-0.0841 

-2.53E-03 

-104 

0 

1 

F 

212 

0.0003 

118 

0.0002 

-0.172 

-1 .00E-03 

63 

39 

2 

B 

237 

0.0106 

121 

0.0076 

-0.174 

-1 .25E-03 

77 

47 

2 

C 

-350 

-0.0106 

17 

0.0076 

-0.267 

-2.32E-03 

-113 

0 

2 

D 

270 

0.0106 

-147 

-0.0076 

-0.164 

-5.39E-02 

114 

-81 

2 

E 

-370 

-0.0106 

-11 




-168 

0 

2 

F 

262 

0.0003 

158 




127 

90 

3 

B 

284 

0.0106 1 

146 

0.0076 

-0.12 

-1.04E-03 

114 

100 

3 

C 

-396 

-0.0106 

17 

0.0076 

-0.214 

-9.08E-03 

-160 

0 

3 

D 

290 

0.0106 




-2.20E-03 

111 

-73 

3 

E 

-407 

-0.0106 

-25 



-1.83E-03 

-120 

0 

3 

F 

264 

0.0003 


0.0002 

-0.246 

-1.41E-03 

104 

89 


Table B.7: 316 Stainless Steel, Path EQ, 610-05. 
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Cycle 


a (MPa) 

e (m/m) 

a 

d 

all 

a12 

1 


0 

0 

1.12 

0 

-146 

6 

1 

c 

796 

0.015 

0.709 

-1.33E-04 

221 

0 

1 

E 

KB 

-0.0002 

0.666 

-6.60E-03 

-255 

25 

1 

G 

-871 

-0.015 


-7.00E-03 

-380 

-25 

1 

1 

760 

0.002 


-1.28E-03 

110 

0 

2 

C 




-1.56E-03 

144 

0 

2 

E 

-800 

-0.0002 

0.152 

-2.70E-01 

-350 

0 

2 

G 

-880 

-0.015 

0.176 

-1.56E-01 

-376 

15 

2 

1 

756 

0.002 

0.25 

-3.12E-04 

135 

0 


Table B.8: Inconel 718, Path I, IN- 16. 


Cycle 

Point 

x (MPa) 

Y (m/m) 

a 

d 

all 

a12 

1 

A 

0 

0 

0.8914 

0 

-145 

0 

1 

C 

535 

1 .06E-02 

0.4399 

-1.59E-03 

-145 

160 

1 

E 

-470 

-1 .70E-03 

0.427 

-1.60E-03 

-145 

-100 

1 

G 


-1 .06E-02 





1 

1 


1 .70E-03 

0.369 

-5.55E-03 

-145 


2 

C 


1 .06E-02 

0.2494 

-3.60E-03 

-145 


2 

E 


-1 .70E-03 

0.303 

-1.59E-03 

-145 

-63 

2 

G 


-1 .06E-02 

0.249 

-2.30E-03 

-145 

-110 

2 

1 

466 

1 .70E-03 

0.302 

-3.08E-03 

-145 

80 


Table B.9: Inconel 718, Path H, IN-23. 
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Cycle 

Point 

o (MPa) 

e (m/m) 

t (MPa) 

y (m/m) 

a 

d 

all 

a12 

1 

A 

0 

0 

0 

0 

0.4032 

0 

-160 

0 

1 

B 

760 

0.0106 

385 

0.0077 

0.138 

-2.22E-03 

144 

85 

1 

C 

-1050 

-0.0106 

127 



-0.99 

-418 

80 

1 

D 

736 

0.0106 

-470 

-0.0077 

0.0959 

-1.45E-03 

-165 

-180 


Table B.10: Inconel 718, Path m, IN-27. 


Cycle 




X (MPa) 

y (m/m) 

a 

d 

all 

mm 

1 

A 

0 

0 

0 

0 

0.999 

0 


0 

1 

C 

740 

0.02 

245 

0.015 

0.1732 

-4.19E-03 

100 

38 

1 

D 

-585 

0.0058 

-400 

0 

0.1239 

-2.07E-03 

-160 

-129 

1 

E 

-857 

-0.002 

-128 

0 

0.186 

-8.49E-03 

-260 

10 


Table B.11: Inconel 718, Path IV, EM-1 1. 
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